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STRESSES  AND  DISPLACEMENTS  IN  A  FOUR  LAYERED 
SYSTEM  WITH  FIXED  BOTTOM 


\  Intermediate  report. 

Contract:  DAJA-85-C-0013  April  30,  1985. 

J 

Introduction.  ^ac^t 

The ‘research,  work  necessary  to  fullflll  the  requirements  of  the  contract 
that  must  lead  to  the  establishment  of  a  computer  program  able  to 
calculate  all  stresses  and  displacements  In  a  four  layered  system  with 
fixed  bottom  submitted  to  a  series  of  loads ,  Is  based  on: 

-  existing  material:  isotropic  multilayer  theory  (BURMISTER,  1943) 
and  anisotropic  multilayer  theory  (y^N  CAUWELAERT,  1983); 

-  original  research  work:  interface  conditions  (fixed  bottom,  partial 
friction)  and  sati sf actor^etfnvergency  at  the  surface  and  in  the 
first  layer  of  the  system. 

This  intermedlat^rgport  will  only  deal  with  the  original  research  work 
that  had^t^HJeperformed .  The  required  original  research  is  completely 
terminated,  which  jusitifies  this  report,  and  at  an  entirely  satisfactory 
level  as  will  be  shown. 

This  report  contains  three  parts: 


complete  discussion  of  the  interface  conditions. 

j 

^/£)The  demonstration  that  satisfying  convergency  can  be  obtained,  and  in 
the  meantime  overflow  problems  can  be  eliminated,  if  the  equations  are 
written  in  closeform  alljipugh  in  a  sufficient  comprehensive  form  so  that 
the  whole  problem  can  still  be  overlooked  but  In  such  a  way  that  all 
numerical  problems  can  be  solved^. 

a^jThe  complete  algebraical  analysis  leading  to  the  equations  In  extended 
form  for  an  isotropic  four-layered  system.  The  analysis  is  also  developed 
for  a  three-layered  system  to  enable  us  to  compare  and  to  check  the 
results  with  those  obtained  by  means  of  reliable  existing  programs.  ^ 


1.  The  interface  conditions. 


The  stresses  and  displacements  in  a  layer  of  a  multilayered  system  are  obtained 
for  an  isotropic  body  from  following  stress  function  (BURMISTER,  1943): 

<j>-  J.C 


and  are_  given  by 


_  C,^  (a-  2y  t  -  TV* 

=  -  jpa  j  3^(jrv*r).  T,  (*v»  ^A\ro  e-  +  ® 

+  C,  w  (a*  1\^\  "»*)  CM+V  « _mz)  e 


1  -w»z 


6  -  WlZ.  Tx-wt^-^Z^t  "|  elk*. 

♦  Ci»v.C>H-wix)e  -3>.v"(^  '  J 


X„ .  -^rVx-yit-)  i>~—  *~v~‘ 

w„  j [A^“  - 

u_ .  .-likibo.  3;cwO  .T,(.^-)  r  ,.  W 

Ei  rvi  L 

-  3>;~  (--»>*■) £"M  <*»" 


I.R.3 


where 

a.  is  the  radius  of  a  uniformly  distributed  circular  load 

p  is  the  value  of  the  vertical  pressure 

r  is  the  horizontal  distance  from  the  axle  in  a  cylindrical  coordinate  system 

Z  is  the  depth 

(J^  is  the  vertical  stress 

is  the  horizontal  radial  stress 
Tr*  is  the  shearstress 

W  is  the  vertical  deflection 

U.  is  the  radial  (horizontal)  displacement 

E;  is  the  Young  modulus  of  the  concerned  layer 
^a.;  is  Poisson's  ratio, of  the  concerned  layer 

are  unknown  parameters  to  be  determined  by  the  boundary  conditions 
Jo  is  the  Bessel  function  of  the  first  kind  of  order  zero 

T,  is  the  Bessel  function  of  the  first  kind  of  order  one 

rr»  is  an  Integrating  parameter 


In  the  case  on  an  anisotropic  body  they  are  obtained  from  (VAN  CAUWELAERT,  1983): 


This  stressfunction  differs  fundamentally  from  the  preceeding  one:  indeed  in 
putting  s.  -  1  in  it,  we  do  not  obtain  the  stress  function  for  the  isotropic  case. 
We  conclude  that  the  two  cases  must  be  handled  in  a  separate  way. 

The  stresses  and  displacements  are  given  by 
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is  the  degree  of  anisotropy,  the  ratio  between  the  vertical 
and  the  horizontal  Young  modulus  of  the  concerned  layer 
is  Poisson's  ratio  expressing  a  strain  in  the  horizontal  plane 
induced  by  a  stress  in  the  vertical  direction 


is  the  index  of  anisotropy. 

1.1.  The  partial  friction  condition. 

Let  us  consider  a  n-layered  system,  consisting  in  (n  -  1)  layers  of  a  finite 
thickness  built  on  a  semi-infinite  body. 

For  each  layer  exists  a  stress  function  fee*6  ,-CjPi^with  4  unknown  parameters: 
the  total  of  unknown  parameters  is  4n. 

Two  parameters  depend  on  the  shape  of  the  load  at  the  surface 

<Tx  =  r^CK 

Trx  *  ° 


where 


«-  c%u 


-m 


At  infinite  depth  stresses  and  displacements  must  vanish  and  thus  An  and  Cn-°- 
We  remain  with  4n  -  4  =  4(n  -  1)  parameters  to  be  determined  with  4  conditions 
at  each  interface. 

The  hypothesis  is  introduced  at  this  stage  that  under  effect  of  the  load,  the 
layers  remain  individually  fully  in  contact,  which  is  expressed  by  imponing  taht 
at  the  bottom  of  each  layer  and  at  the  surface  of  next  layer  vertical  stresses 
(<TX),  shearstresses  ("Cr*  )  and  vertical  displacements  (W)  are  identic. 

The  fourth  interface  condition  depends  on  the  relative  adhesion  in  the  horizontal 
plane  between  the  considered  layers. 

The  two  extremes  are 

-  full  continuity,  expressed  by  setting  that  the  horizontal  displacements  (u) 
are  identic; 
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-  frictionless  interface,  by  considering  the  interface  as  a  principal  plane 
and  thus  by  setting  the  shearstresses  equal  zero. 

Partial  adhesion  has  been  temptatively  introduced  by  several  authors,  utilizing, 
in  the  same  way  as  WESTERGAARD  (1926),  a  relation  between  horizontal  displace¬ 
ments  and  shearstress: 


K( 


lc\ 


U.  \ 


)-  XrZi 


where  u^  is  the  horizontal  displacement  at  the  bottom  of  the  i-th  layer  and  u.j+1 
that  at  the  surface  of  the  (i  +  l)-th  layer. 

We  shall  prove  that  such  a  relation  cannot  be  correct  in  the  case  of  a  multilayer. 

One  has,  for  an  isotropic  body,  following  relations  between  displacements, 
shearstrains  and  shearstresses: 

aw 
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We  know  from  the  boundary  conditions  that 

Wi  «  »  TrZ-.-n 

This  is  true  evetjwhere  on  the  interface  so  that  we  also  can  write  that 

r&r  ~ 

By  substraction  we  obtain 
fbu.-,  _  ^;±!_  ~ 
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We  can  thus  write 
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This  relation  must  be  satisfied  for  all  values  of  the  parameter  m,  so  that 
one  must  necessarily  have  that 
.  u_i 


"te7 


*1  f&T.  S' 

The  solutions  of  those  differential  equations  are 

zVA 
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Comparing  those  solutions  with  the  relation  above  for  the  horizontal  displacements 
we  conclude  that  the  obtained  expressions  must  be  deduced  from  a  stressfunction 
different  from  the  original  one  which  is  nevertheless  the  unique  solution  of  the 
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compatibility  equations.  Comp  atibility  is  thus  not  respected  and  relation 
■  trz  cannot  be  accepted. 


Our  meaning  is  that  the  only  way  to  express  partial  continuity  (or  partial 
adhesion)  consists  in  writing 


with  fe  £• 

When  k  =  1,  one  has  full  continuity 
When  k  t  1,  one  has  partial  continuity. 

It  is  necessery  now  to  give  a  physical  sence  to  the  parameter  k. 

Excepted  the  extreme  case  of  a  frictionless  interface,  there  will  always  be 
some  friction  between  the  layers  at  their  interface. 

We  rely  then  our  approach  on  Coulomb's  definition  of  friction.  If  is  the 
angle  of  friction  at  the  interface,  there  will  be  no  sliding  (in  the  geotechnical 
sense  of  the  word)  between  the  two  layers  as  long  as  Trz-<.<3'*  M 
The  limit  value  for  k,  at  which  sliding  due  to  shearfailure  will  occur,  is  then 
given  by  -  |v^>. 

Beyond  this  value  of  k,  shearstresses  vanish  and  the  interface  has  to  be 
considered  as  frictionless. 

We  must  of  course  take  into  account  here  the  stresses  due  to  the  wheel  load 
but  also  those  due  to  the  own  weight  of  the  layers  above  the  interface,  which 
reduce  to  vertical  stresses  only. 

The  final  relation  becomes  then 


‘Cfa  ^'rr<x\^c) 


The  values  of  the  stresses  due  to  traffic  varie  with  the  distance  to  the  axle 
of  the  load. 

We  carry  then  the  calculations  out  in  two  steps: 

-  we  calculate  for  different  values  of  k  the  maximum  value  (function  from 
the  distance  from  the  axle  of  the  load)  of  the  ratio  T/T'. 

-  we  determine  the  limit  value  of  k  in  function  of  tgi^. 

As  an  illustration  of  the  method,  we  have  taken  the  most  simple  case,  that  of  an 
isotropic  two-layer  system,  with  H  the  thickness  of  the  first  layer. 

For  simplicity  we  take  u.,  »  =  o.5 
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We  write 
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The  boundary  conditions  are  then: 

At  the  surface  (z  *  0): 

^2.-=  ^ 
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At  the  interface  (z  =  H); 
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Solving  the  system  of  6  equations,  one  obtains  at  the  interface  (z  =  H): 

r  <0 

tT2*  VpT.O-^T.C^O  [t*-*0^",rt  4  (-.♦?)£*""] 

V  a  ^-F)^-L)e^H4  l^A-TX.)  -a^H(f-L)  +  ^H,(xi-Fl)yi-H 

t  ^>I+F  )  ^4  L) 


I.R.8 


We  have  performed  the  computations  for  different  values  of  H/a,  with  a  =  10  cm, 
2 

p  =  0,6  MN/m  and  *  10.  We  have  considered  a  specific  weight  of  the 

first  layer  of  22  kN/rn*. 

The  results  ar  given  on  figure  1:  in  abscissa  one  finds  the  values  of  the 
ra tio’Vtf'"  and  in  ordinate  the  relative  thicknesses. 

The  curves  give  for  different  values  of  k  the  maximum  value  of  Ty/<p . 


I.R.22 


i  i 


This  is  practically  impossible. 

One  could  of  course  be  tempted  to  interrupt  the  integration  procedure  when  the 
first  integral  has  converged  to  a  satisfactory  level,  "hoping"  that  the  second 
integral  can  be  neglected  at  that  moment. 

To  illustrate  the  danger  of  such  an  approach  we  return  to  the  semi-infinite  body. 

In  the  case  on  an  isotropic  body,  the  deflection  at  the  surface  and  the  vertical 
stress  at  a  depth  z  are  given  in  the  axle  of  the  load  by 
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Those  integrals  can  of  course  be  solved  analytically 
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We  can  compare  (2)  with  the  second  integral  of  (1)  and  (3)  with  the  first 
integral  of  (1) . 

We  perform  then  a  numerical  integration  of  (2)  and  (3)  and  stop  the  procedure 
when  (3)  has  converged  to  a  satisfactory  level,  which  is  easily  checked  by 
comparing  the  obtained  result  with  the  correct  one  given  by  (5). 

The  difference  between  the  numerical  result  for  w  obtained  at  that  moment 
by  integration  of  (2)  with  the  analytical  result  given  by  (4)  will  give  us  an 
illustration  of  the  possible  error  when  integrating  (1)  and  stopping  the 
process  when  its  first  integral  has  converged. 

This  difference  is  illustrated  in  figure  2. 

In  abciss,  we  have  the  convergence  level  adopted  for  the  vertical  stress  and  in 
ordinate  the  error,  expressed  in  %,  on  the  values  of  the  vertical  stress  and 
the  vertical  deflection. 

_3 

One  sees  that  for  even  such  low  levels  as  10  ,  the  value  of  the  vertical  stress 
is  absolutely  correct,  while  the  error  on  the  deflection  varies  between  -  5% 
and  +  8%,  depending  on  thechosen  convergency  level  and  the  relative  depth  at 
which  the  vertical  stress  is  computed;  worst  of  all  is  that  we  have  no 
means  to  predict  either  the  direction  either  the  amplitude  of  the  error  on  w. 


\ 


2.3.  The  vertical  deflection  at  the  surface. 


When  one  computes  the  deflections  at  the  surface,  convergency  is  obtained 
only  very  slowly. 


only  very  slowly. 

To  illustrate  this  let  us  look  at  the  expression  of  the  deflecti 
surface  developed  in  §  2.1. 
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For  large  values  of  m,  numerator  and  denominator  tend  both  to  F,  so  that  for, 
let  us  say  m  =  m^,  the  expression  above  could  be  written  as  follows: 
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The  first  integral  converges  fast,  the  second  converges  proportionnaly  to  1/m. 
This  means  that  if  one  should  want  a  result  correct  at  10-5,  one  has  to  perform 
the  numerical  integration  of  the  second  integral  until  values  of  m  above  100,000 
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Underflow  will  obviously  occur  now,  but  most  of  the  computers  have  a  routine 
that  sets  variables  subjected  to  underflow  equal  to  zero.  If  such  a  routine 
does  not  exist,  it  is  very  easy  to  build  it  into  the  program. 

But  more  interesting  is  the  fact  that,  having  transformed  the  relations  for 
Cj  and  D^,  convergency  will  occure  quite  quickly  and  in  a  complete  safe  way: 
the  numerators  both  tend  to  zero,  while  the  denominator  tends  to  a  constant  F 


This  can  be  obtained  automatically  in  writing  the  boundary  counditions  at  the 
surface  (z  =  -H)  as  follows: 


.  -SrviH  -h*H 

A,e_  *  b.e,  - 


However  this  is  only  true  in  the  case  of  a  two-layer  system. 

In  a  three-layer  with  H^,  the  thickness  of  the  first  layer,  and  H^,  the 
thickness  of  the  second  layer,  occur  exponents  such  as 

But  they  eliminate  when  writing  the  denominator  in  closeform  so  that  dividing 

nmU  o_u  * 

the  expressions  by  the  largest  out  of  e  l  and  e  2  is  e  nough.  If  one 
should  divide  by  e2mHi*2mH;>  >  the  denominator  would  also  tend  to  zero,  which 
should  stop  the  program  because  of  dividing  by  zero. 

Here  is  another  reason  for  writing  the  equations  in  closeform  for  three- 
and  more- layers. 
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2.2.  Over  and  underflow  problems. 

Dur.ing  the  integration  procedure  m  varies  from  0  to  a  value  high  enough  to 
ensure  convergency.  We  mean  by  this  that  the  integration  procedure  can  be 
stopped  from  the  moment  on  that  the  terms  of  the  series  become  so  small  that 
they  have  no  more  influence  on  the  final  result  and  can  thus  be  neglected. 

Practically,  however  this  means  that  m  can  reach  quite  high  values  such  as 
20  or  30  for  example. 

To  illustrate  the  influence  of  this,  let  us  go  back  to  the  two-layer 
developed  in  the  preceeding  paragraph. 

The  values  of  and  D^,  from  which  the  values  of  all  the  other  parameters 
can  be  deduced,  are 


The  geometrical  unities  are  generally  expressed  in  function  of  a,  the  radius 
of  the  load. 

Let  us  consider  H/a  =  5. 

One  immediately  sees  that  no  computer  can  handle  exponents  as  emH//a  and  e2mH//a 
without  overflow  occuring  for  values  of  m  above  10. 

However  this  problem  can  easily  be  solved  by  dividing  both  numerator  and 
denominator  by  e  : 


C»  * 


-  F  vtvi  n  Je.  -  £  A  -  i  )  c. 


H] 


F  *  -  (»+  Vh1)] 


H  _\  -  f) 


A 


_  (F- *-») 
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Replacing  Aj,  B^,  Cj  and  Dj  by  their  values,  the  deflection  becomes 

w. "iMlM 

e.  1.  *> 


[i—  H 

«=*77 


"  QF-^- 

(2.F  -*  2w>  H 


awt4)  1 


At  the  origin  of  the  integration  (m  =  0),  the  term  in  between  brackets 
becomes  indefinite:  0/0. 

This  has  no  influence  when  computing  stresses,  because  the  Bessel  functions 
products  occuring  here  are  also  zero  at  the  origin:  J0(mr)  .O^ma)  =  0 
for  m  =  0. 

But  in  the  case  of  the  deflection 

K  l>C™r)  .  _  ck 

It  is  therefore  absolutely  necessary  to  have  the  term  in  between  brackets 
in  close  form  to  be  able  to  determine  its  value  for  m  =  0.. 

The  importance  of  the  first  term  of  the  series  is  not  negligible:  for  m  =  0 
the  term  in  between  brackets  is  equal  to  E*  /  [£j.  . 

If  h  is  the  interval  choosen  for  the  numerical  integration,  one  can  then  write 


—  J  _L. 

»  |  z  3  1  J  ) 


and,  if  we  make  a  semi-infinite  body  from  the  two-layer  (E1  =  E,,,  ) : 

W j  T  *  f  3^1-)  L 

Comparing  this  expression  with  that  for  the  deflection  at  the  surface  of 
a  semi -infinite  body 

W«o  =  -  /£ 

one  concludes  that  the  contribution  of  the  first  term  h/6  is  indeed  not  negli¬ 
gible,  especially  when  we  have  in  mind  that  the  only  practical  measurement 
that  can  be  performed  on  a  real  roadstructure  is  the  vertical  deflection 
at  the  surface. 
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At  the  surface  (z  =  -H): 

A,e"'H+  D,  o- 

A,el~M-  0,«rH  *  t>,  (^..wH)^^ 

At  the  frictionless  interface  (z  =  0): 

A,+  1b, -C,  (4-2^0+  ,  &,.+  Da 

A,~  i*>,  4  1) ia,CL,  V  ,  sr  «o 

~  +■  2fn.  T»-3_=  o 

-  i<=,  OVJ-iM-'-Vi] 

=  [-  &2  _  1T>J  V)J 

Solving  the  system  for  and  0^,  one  obtains  (BURMISTER,  1943) 

C,=  ^(/I-  F  .wH)«TH  _  -^- 

(f-wH>TM].  CL_ 

where 

V-  +  (zF--0.2h*H 


F-  (*-)-0 


Er  [»•*•>- 0 
E*  (x+>av) 


A,  •  c,  o  *K>  -  ■», 

B(  „  C,F  -  ‘D,['l-l+''-'F) 


The  vertical  deflection  at  the  surface  is 


a  +  Ui 

!>*•  e7 


r«fl"Q  T,fa.«)  J  ft, we 

Ttr,  L- 

.  (1-4)*,-  v-  H)  T>,u^H]  eU* 
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2.  Solution  of  particular  numerical  problems  (convergency-preblems) . 

2.1  The  full  slip  interface  condition. 

The  value  of  any  stress  or  displacement  is  obtained  from  one  of  the  above 
mentioned  relations. 

Let  us  consider,  for  example,  the  vertical  stress  in  the  i-th  layer  of 
an  anisotropic  layer: 

ffj  »  ^pcx.  ^  0r*r')  •  C'vl  *■)  (*■*■  £A|Vv*  e-  +  3 

•V  ^CiSi  W»V«f*VVl  4- 

The  integration  can  only  be  performed  numerically. 

Thus  one  must  calculate  the  value  of  the  stress  for  a  set  of  values  of  m 
growing  from  0  to  a  value  high  enough  to  ensure  converaency. 

For  each  value  of  m,  those  of  the  parameters  Ai ,  B.. ,  and  must  be 
determinated  out  of  the  set  of  boundary  counditions,  a  system  of  (4n  -  1) 
equations  with  (4n  -  1)  unknowns  in  the  case  of  a  fixed  bottom  and  n  layers 
above  it. 

The  first  programs  solved  this  problem  by  inverting  the  matrix  of  the 
(4n  -  1)  unknowns.  Nevertheless  the  inversion  procedure  leads  in  some 
cases  to  unsoluble  difficulties  because  of  the  presence  of  the  negative 
exponents  tending  to  zero  in  the  determinant  of  the  denominator. 

Other  programs  have  tried  to  avoid  the  inversionprocedure  as  follows:  one 
chooses  appropriate  values  for  Bn  and  Dn>  goes  trough  the  whole  set  of 
equations  and  verifies  in  how  far  the  surface  conditions  are  met.  One 
then  chooses  another  pair  of  values  for  Bp  and  Dn  and  follows  the  same 
procedure.  Since  the  whole  process  is  linear,  the  correct  values  for  Bp  and 
Dn  can  finally  be  obtained  by  linear  interpolation  after  two  runs.  The 
difficulty  lies  in  the  appropriate  choice  of  the  values  of  Bp  ard  Pn  to 
ensure  a  numerically  correct  interpolation. 

However,  even  those  programs  are  not  entirely  appropriate  for  the  cases  with 
frictionless  conditions  at  some  interfaces. 

We  shall  show  this  with  the  most  simple  case,  that  of  a  two  layer. 

Writing  A,,t>j,C{  ,D  {  instead  of  A,vv»\ ^ivvi* the  boundary 
conditions  are  in  the  case  of  two  isotropic  layers,  with  the  origin  (z  =  0) 
at  the  interface,  the  thickness  of  the  first  layer  being  H  and  the  second 
layer  semi -infinite: 
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We  also  conclude  that  the  relative  Influence  on  the  deflection  Is  much  more 
Important  when  we  fix  the  horizontal  displacements  (u  =  0).  This  should 
be  the  case  In  a  laboratory  testplt  with  lateral  walls,  but  less  In  the 
case  of  a  real  road  where  lateral  movements  are  not  restricted. 

The  relative  Influence  of  the  condition ta=  0  Is  also  more  Important  than 
that  of  the  condition  w  =  0,  altough  less  Important  than  the  condition  u  =  0. 
It  seems  nevertheless  very  unlikely  that  there  would  be  no  friction  between 
the  subground  and  the  last  layer. 

The  easiest  way  to  fix  the  bottom  from  a  mathematical  point  of  view  is  on 
the  other  hand  the  condition  w  =  0. 

Taking  then  Into  account  the  little  influence  of  the  chosen  condition  on  the 
deflection  at  the  surface  and  the  fact  that  conditions  u  =  0  and  trt=  0 
have  less  physical  sense,  we  shall  retain  the  condition  w  =  0  as  the  most 
indicated  fixed  bottom  condition. 
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The  deflection  at  the  surface  is  given  by  the  same  relation  as  above  with  the 
appropriate  value  for  C. 


One  sees  that  in  the  3  cases,  the  deflection  is  composed  of  a  first  term 


r 


E  b~S)  g  ^ 


cLwy 


.  sf-A-vy) 


This  term  is  the  deflection  on  top  of  a  semi -infinite  body. 
In  the  case  of  an  isotropic  body,  one  has  (n  =  1) : 


The  second  term  wr  depends  on  the  choosen  boundary  condition;  but  In 
the  three  cases  it  reduces  the  value  of  because  of  the  fixed  bottom* 


We  have  computed  the  values  of  and  for  different  values  of  H/a.  The 
results  arg^  given  below,  at  a  factor .  for  s  =  0.5  and  for  the  isotropic  case 


s  =  0.5 

w  =  1.025 

H/a 

wr  (T  rz  -  °) 

wr  (■**  °) 

wr  Os 

1 

0.339 

0.729 

0.153 

2 

0.120 

0.295 

0.043 

3 

0.058 

0.147 

0.019 

4 

0.034 

0.086 

0.011 

5 

0.022 

0.057 

0.007 

6 

0.015 

0.040 

0.005 

7 

0.011 

0.029 

0.004 

8 

0.009 

0.023 

0.003 

9 

0.007 

0.018 

0.002 

10 

0.006 

0.015 

0.002 

n  =  1 

=  1 . 000 

1 

0.402 

0.934 

0.276 

2 

0.143 

0.442 

0.086 

3 

0.069 

0.231 

0.040 

4 

0.040 

0.138 

0.023 

5 

0.026 

0.091 

0.015 

6 

0.018 

0.064 

0.010 

7 

0.013 

0.048 

0.008 

8 

0.010 

0.037 

0.006 

9 

0.008 

0.029 

0.005 

10 

0.007 

0.024 

0.004 

We  conclude  that  from  a  depth  of  about  H/a  =  5,  the  absolute  Influence  of  the 
fixed  boundary  is  negligible.  This  influence  will  still  be  much  lesser  in  the 

case  of  a  roadstructure  where  the  E-modull  of  the  layers  are  sensitively  higher 

than  the  modulus  of  the  subground. 
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We  now  choose  an  appropriate  function  so  that  w  becomes  zero  for  z  =  H. 

\(r)  -  -  )>«■  I  -  e>r»V“h) 

O 

l*+Y)  1  J 

The  final  expression  for  the  deflection  is  then 

W  =  W,  4  \{r) 

One  verifies  that  is  indeed  only  a  function  of  r  and  that  w  =  0  for  z  -  H. 

For  the  same  reasons  as  those  developed  in  §  1.1.2,  one  of  the  parameters 
A  or  C  must  be  zero. 

The  other  paranster  is  obtained  by  a  supplementary  boundary  condition 
Ca  mechanical  condition): 

-  TV*  *  o  at  the  depth  H 

U.  =o  at  the  depth  H. 

If  we  still  suppose  and  thus  A  =  0,  one  obtains  in  the  case  that  Trz  =  ° 

ns(w^)  Cf  ha2  »  - - - — - — - 

-  (*-s)  - 

The  deflection  at  the  surface  is  then 


.sM-*) 

VV»  “IT7 - \~ 

F  E  i^-s) 


^  i**"-  ”c* 


_w,W 


r  r-  z\  2s(v»-0  .  JL* — e 

♦  (*+)<*)  Wrv,  S 


In  the  case  u  «  0: 


e',Mm 

ns(m>~)  Cs~'  ji7(W>« 


-2  .fW*, 
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1  ^2_.  2^_Fi  xed_bottom_exgressed_b^^!pe£t’§!3iC33  _COPdt  tiQD_QQly 


Referring  to  the  stress  and  displacements  equations  given  in  paragraph  1,  the 
condition  w  =  0  at  a  depth  H  is  written 

V  r^^\cs^sm  -  WeJ>"  )  =  0 


Replacing  B  and  D  by  their  values 


i  if  ,  ^  -V"»  .  Zr>s(y>^y)  a~s^‘\ 

>  -J^TT  J 


4Dw,LW  '  j 


■S^Vj  _  A —  e 

(*»V0 


-wH 


During  the  integration  proces,  the  value  of  m  tends  to  infinity,  so  that  the 
limit  expression  of  the  equation  becomes 

km  Y>Ul  Crw'sC^H1®^] 

*0  *“ 

This  relation  can  only  be  satisfied  by  setting 


-ajw'i 


Cl  -  o  when  ■*  >  'l 

^  —  o  when  -s  4 

It  is  easily  shown  that  in  the  case  of  an  isotropic  body,  one  must  always  have 
C  =  0,  because  of  the  factor  z  multiplying  C  in  the  stressf unction. 

Taking  s<4  ,  C  becomes  _  H  . 

e‘mH[0+We  J 

r>s(*+y)  Csw1  — -  u  _JW,U  ~  "  '  w  ..  v  -?smn 

v-  2  Cv,4-H e 

The  expression  of  the  deflection  at  the  surface  and  in  the  axle  of  the 
load[3%0,v,r) for  r=«J  is 

li?;?i.El?§d_hottom_by_a-_2eometri  cnl_and  _a_me  chan  i  cal  _conditioni 
The  vertical  deflection  at  a  depth  z  is  given  by 

w, .  t>>~v  “••>-) 

Vjsfvt^)1  fc swiV1*11-  3>$vwVJI*,x)"\ 
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But  with  the  general  solutions  of  the  compatibility  equations  In  multilayer 
theory,  we  can  do  It  also  In  another  way  by  expressing  that  w  *>  0  at  the 
desired  depth  and  determinating  the  corresponding  values  of  the  parameters 

V  Bn*  Cn  and  V 

Since  there  are  thus  several  possibilities  to  express  a  same  boundary 
condition,  It  is  necessary  to  compare  the  results  obtained  and  retain  the 
one  that  seems  the  most  appropriate. 

To  do  this  we  shall  consider  the  most  simple  case,  that  of  the  semi-infinite 
body:  the  one-layer  case. 

1 j2il^_Baslc_eguatIonsi 

Let  us  consider  a  semi-infinite  anisotropic  body  submitted  to  a  uniformly 
distributed  vertical  pressure  at  its  surface. 

The  stressfunction  is 

The  surface  boundary  conditions  >T«x-so  ,  z«ro  )  are  deduced  from 

the  relations  given  in  §  1.  for  the  stresses 

Bm*1)  ♦  I’stvi  }-  -i. 

Solving  this  system  for  B  and  0,  one  obtains 

IWs  -  x  +  f  2wj£*n>*)Cx*M 

r\  (A-s)(*+y)  D**?  =  -  2«  n 

The  next  step  depends  on  the  boundary  condition  that  fixes  the  deflection 
at  a  depth  H. 
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1.2.  The  fixed  bottom  cSndition. 

The  boudary  conditions  discussed  in  previous  allnea  impHcate  that  the  last 
layer  of  the  multilayer  is  considered  as  a  semi-infinite  body. 

One  can  also  consider  the  case  of  a  multilayer  built  on  an  undeformable  body, 
that  thus  any  vertical  displacement  vanishes  at  the  contact  f®ce  w'th  the 
undeformable  body:  we  shall  call  this  a  fixed  bottom  condition. 

This  condition  can  be  introduced  in  several  manners  and  thus  demands  a 
detailed  analysis. 

A  vertical  displacement  is  obtained  by  integration  of  the  vertical  strain: 

W  r  J  t*  ■<** 

It  would  not  be  correct  to  resort  to  an  integration  between  limits,  such  as 


where  H  could,  for  example,  be  the  depth  at  which  we  want  the  bottom  to  be  fixed. 
In  doing  so.  we  would  ignore  the  contribution  (zero  or  not)  to  the  vertical 
deflection  (or  displacement)  due  to  other  parts  of  the  body  that  we  neglect 
by  integrating  between  specified  limits. 

The  correct  way  consists  in  writing  (TIMOSHENKO,  1970): 

W= 

where  is  a  fonction  of  r  only,  and  thus  a  constant  regarding  Z  ,  so  that 
by  differentiating  we  obtain  again 

^  -  Ez 

By  choosing  an  appropriate  expression  for ^(r)  ,  we  can  obtain  the  bottom  fixed 
at  the  desired  depth;  by  doing  this,  we  introduce  in  fact  a  geometrical 
condition  fixing  the  reference  level  for  the  vertical  deflections  at  the 
choosen  depth. 


Let  us  consider  an  average  value  of  tg'p  =  1.5,  value  utilized  in  the  design  of 
continuously  reinforced  concrete  pavements  (Me  CULLOUGH,  1981). 

We  deduce  from  figure  1  the  limit  values  for  k: 


H/a  =  1  k  =  0.20 
H/a  =  2  k  =  0.17 
H/a  =  3  k  =  0.10 
H/a  =  4  k  =  0.03 


If  we  consider  a  surface  layer  built 
the  limit  values  become: 


on  a  sand  basecourse 


(\^=3o°,  a.5%  ) , 


H/a  =  1  k  =  1 

H/a  *  2  k  =  0.60 

H/a  =  3  k  =  0.40 

H/a  =  4  k  =  0.25 


With  our  approach,  the  case  of  full  continuity  (ui  =  ui+1,  often  called  full 
friction)  becomes  a  particular  case  for  which  the  angle  of  friction  between  the 
layers  corresponds  with  a  value  of  k  =  1.  For  other  values  of  k  we  have 
partial  continuity  (u^  /  u-j+i»  what  c°uld  he  called  partial  friction  or  partial 
adhesion) . 

We  notice  that  for  values  of  k  smaller  than  one,  u..  is  smaller  than  u..+j. 

This  means  practically  that  the  lateral  movements  of  the  surface  layer  are 
retained,  for  example  by  shoulders.  Such  a  construction  reduces  the  vertical 
stresses  on  the  subground  which  improves  the  lifetime  of  the  road  structure. 

For  values  of  k  larger  than  one,  u^  is  also  larger  than  u^+^.  Here  the 
lateral  movements  of  the  surface  layer  are  easier  than  those  of  the  sublayer, 
as  in  overlayconstructions  for  example. 

For  the  limit  value  k  =  ,  one  obtains  alsoTf*=  0. 

This  case  is  called  the  frictionless  interface. 


For  k  =  *0  ,  one  obtains  indeed 


_£(  -2v^vi’E\] 


-  •  cIk* 


which  is  the  value  for  in  the  case  of  a  frictionless  interface  (BURMISTER,  1943) 
as  will  be  shown  in  §  2.1. 


Fiqure 


I.R.24. 


The  only  way  to  solve  the  problem  satisfactory  is  to  split  the  expression 
of  the  deflection  in  another  way  than  the  one  we  had  done. 

We  first  write  the  expression  of  the  deflection  with  negative  exponents  only: 


W  a  -}}C\. 


E, 


9k 


X,  (mr) .  £**>») 

Yvn 


l  F  +  [C2F-^2w,H  -  *  (>i-F)eWl  j 


We  then  divide  the  numerator  of  the  term  in  between  brackets  by  the  denominator: 

2 (a-  |  X(wtr)  .  3,(w<i  fc) 


W=  -|i< 


r* 


p+  [(2f-))2w,h  -  -*•  n  t 


A  + 


-4w,» 


and  split  the  integral  into  two  parts  from  which  the  first  is  integrable 
analytically  and  the  second  converges  in  the  usual  way. 


For  r  =  0,  one  has 

Wl 

For  r  =  a,  one  has 
/<* 


.  oLua.  s  A 


I  v-r 

Jo  \m 

For  r<^  a,  one  has 

F(T’--5-i  "*•> •£) 

o 

where  F  is  the  hypergeometric  function  of  GAUSS: 

r*  /  \  ~\  -sr  CtM-OO"  v, 

F  (c,k,c>»)  .  <21 


cA  v»-o 


For  r)  a,  one  has 


C”'1 0  •  T,  [ry 


dlw»- 


The  obtained  result  will  now  be  correct,  while  convergency  is  reached  as 
fast  as  for  the  other  equations  for  stresses. 

But  again,  if  we  are  to  be  able  to  compute  as  indicated,  we  must  have  the 
equations  in  closeform  at  our  disposal,  although  in  such  a  form  that  the 
integral  can  be  split. 


1 
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2.4.  Converqencv  in  the  first  layer. 


As  for  the  deflection  at  the  surface,  the  numerical  computation  of  the 
stresses  in  the  first  layer,  in  fact  nearby  the  surface,  also  converges 
very  slowly. 

To  illustrate  this  let  us  look  at  the  relation  for  the  vertical  stress  in 
the  first  layer  (the  equation  is  given  in  §  1.): 

(Tx=  X (yv'r)  j_Ay*  e.  *  t 

.  c,~  t*-  V. 


We  replace  A^,  ,  C^,  Dj  by  their  values  obtained  in  §  2.1. 


a  (  -W.H  -w«x. 

<T*»  Jscx.  1  F  ( H* *)J 


-w,H  wix 

*-  L  >  v 


L  >  '  ' 


L  3*-  e 


The  values  of  z  are  negative  In  the  first  layer  (z  *  0  at  the  interface). 

The  term 

-  ,  vT  -w*H  -v"'x  A 

converges  very  slowly  for  values  of  -z  nearly  equal  to  H,  the  other  terms 
converge  normally. 

To  solve  the  problem  created  by  the  first  term  we  divide  again  numerator  by 
denominator: 

.  U  >1  r  -VwH  -VviX 


-Ai^M 


The  second  term  of  the  second  member  converges  normally  so  that  we  can  again 
split  the  integral  in  several  parts  from  which  the  one  that  converges  slowly  is 


hev  To0"0- I'i  +'*)]«■ 

r  i  *- 


-VwB  -h«t  , 


/ 
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This  integral  is  known  as  a  LIPSCHITZ-HANKEL  integral,  but  only  some 
particular  cases  are  integrable  analytically.  To  solve  the  problem  for 
all  cases  we  have  to  make  a  detourthrough  the  analysis  of  stresses  and 
displacements  in  a  semi -infinite  body  submitted  to  an  isolated  local  force  P. 

i _§££§§§§§_ i sglacements_u nder_an_ isolated _ load . 

The  Hankel  transform  in  the  case  of  a  uniformly  distributed  load  is 

3*,  [a-**)  <k  d.Wi 
tv>  ' 

2 

We  consider  the  resulting  load  P  =Tpa  acting  on  a  surface  whose  area 
reduces  to  zero. 


The  relations  for  the  stresses  and  the  displacements  are  then  given,  in  the 
case  of  an  isotropic  body,  by 


cAvvi  -  'Z  Y'n*  ol>»n 
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air  (rW)54 

i=  -  X-  a_hV)  \  %  Ur'W* 


p  af-i-h1)  ^  _  _T__  f»+H  ^ 

a  XTT  E  (aW)'4  air  £  (a\r1)^ 

,  JL  f'-iHMd  l -  £s_  (fid  fL 

aT  £  1  air  E  J, 

_  p  (r\aL)V-2  p  (>»+>-)  *«" 

~  2.T  E  f(rW)Vi  "air  E  '  (*Vr')5A 


2j:4J12i_Stresses_and_di sp2acements_under_a_uni formjy^ di stri buted  load. 

The  stresses  and  displacements  under  a  uniformly  distributed  load  can  be 
obtained  by  integrating  the  relations  under  an  isolated  load  over  the 
concerned  area 


<v- 


air  ,  a. 


fdf 


where  Ox  Is  given  by  one  of  the  relations  of  §  2.4.1.  wherein  the  distance 
r  must  be  replaced  by  £cL4  pl  -  2rj>  us-i©}'4 

?ifi?i.B§l§rions_for_the_stress_0>_ 

The  relation  for  rr  under  a  distributed  load  is  given  by 

fly  =  (nor)  .  7,  (*”  -  ycx  J  Jc(«vtr)  0,  c4ki 

0  .*  *  <<* 

_bA  (>1-^)  ?iW  t,hA  v  ya 

'  Js j  vnr  Aa  r 

so  that  comparing  with  the  expression  for  the  stress  under  an  Isolated  load 
we  can  conclude  that 
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2_.  4i4i_Rel  ati  ons_for_  the_stress_<Ti. 

These  relations  can  be  deduced  from  those  established  for  the  stress  0> 
?i4.;5i_Re1atlon_for_the_stress_T<-z 

Trx  *  7,  (W).  (wv^)  enatVV'XcAwi 

X,..  i  ^  "3J[rVf-2rp^ 

4  27rA  L  (TlVr’v  j>l-2rpun$)  ^ 

?i$i§i_B®l§£i9U§_f2C-£!3§_¥ertical_disglacement  w. 
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aTJo  X  (aVr'vf »7-2rf»wjftV' 
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2. 4. 7i_Rel§ tlons  _for_the_hori zgntal_dlsj2lacement_ui 


I  C-  '  Vo 

*-  >0 

_  Ov  feth?  1  3\  (vor)  .T,  x 


a>  *'  171 1  -  (f  A^-lrj>uv6) Vx  ^>r\^-lr{in  6)  '/u 

Io,:^  ^(v^ryx^^a^dl^r  -t-TP  *  (Vv  f1-  Iff  *»*)' 

'*>  ■'»  ^  ( 2*  +r1  +  f>1  -  2  rf> 


?i4.81_Resolution_of.the_doubleJntegraIs. 


The  integrals  of  alinea  2.4.7.  are  most  easily  solved  in  transforming  the 
variables  &  and  p  by  setting 

^  ~  j>  Y  =  \  hvi 0  At)  =  dx 

One  obtains 

_  *  r  _w^)fc _ j. 


T  .  i_  - — - *  rU 

2F  /  (aVxVr'-axO^Sr^^-ixO^ 

—  Pt 

that  can  easily  be  computed  numerically. 

T  l1  2a  ( rS  y'-'Xrx) ( c\~ x; ) /r  , — 1— - -  *  —  rT"_  > 

->s 

v  v  r  _ * 

2TT  (at+xVr1-lxr')(x2v-f'vc^-2xrV/l 

*  -  W 

-  L  *  Axlc'-^  _ _  * 

**  (o.'.x'+t'-lr*')  (V t 

I3  ar.d  I4  are  particular  lipschitz-Hankel  integrals  that  we  shall  solee  in 
next  alinea. 


I5  can  be  deduced  from  I 2. 
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xt-  JL 

2ir 


(r’-Zxr  +M)*  -  («.■*-*’)  * 


I7  can  be  deduced  from  1^. 

Ig  and  Ig  are  particular  Lipschitz-Hankel  integrals. 

2.4.9.  Resolution  of  the  Lipschitz-Hankel  integrals. 

The  solution  of  the  Lipschitz-Hankel  integrals  is  given  by  WATSON  (1960). 

r*  »■* !,(.»).  V  = 

'o 

_r*lv  V. 1  )  *■  1  1  ' 


TT  oT+  V  r  (2v+  * ) 


where 


to 


T  -  V 

-L-b-  r-p 


V>Vc-  Ik  ^ 


oO 


_  .  ov^  _ 

=  P  1TT- 
with  (WAYLAND.1970) 


IV* 


TT 


T{«.»4s  “•)*] 

and  WS  »*r'-  2.W"'4' 

The  resulting  integral  can  easily  be  solved  numerically. 
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I4  -  T,  (»v>r)  •  T,  cUa*. 

->?4 

*% 

n\^-^  o#rA 

Ig  can  be  deduced  from  Ig 
Ig  can  be  deduced  from  1^; 

?i^iiQi_.l?B!T§§§l9D§_foF_con)gutati9ns_in_the_axle_gf_the_load. 

When  stresses  and  displacements  are  computed  in  the  axle  of  the  load,  one  has 

«•»  - 2. 

bOv-  I  SLY*X  - - ====— 

•'<>  CX.  y*N-x2 

\*  _yV,-X  OT~ 

J-i'  cAa~=  ~^~^pyvT 

x,..  t>-  v>^  r  ^ 
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10*  1-0  y**  TV 
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T  ,  \  X  ( Vw  fc)  ~W!i  . 
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I  °*  ,  N  -W-T  , 

X7.  :  1  X  ( •*»  0  ^  *-  C^Vv'  - 


X^x1  -  Z 


a.  .l-i 


X  2  :  r  =  VO 


Xft  : 


0,  :  )> 


X,  tv^O  T,  (y**)  at  ^cU-  -  o 


To  be  applicable,  all  the  developments  of  paragraph  2.4.  again  require 
all  the  equations  to  be  available  in  closeform. 
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3.  The  complete  algebraical  solution. 

We  shall  develop  here  the  complete  algebraical  analysis  leading  to  the 
different  computerprogramsin which,  because  of  the  particular  method  that 
we  have  adopted,  all  the  problems  detail  led  above  are  solved. 

But  to  make  the  understanding  of  the  method  easier  we  first  develop  the 
analysis,  with  the  necessary  comments,  for  a  three-layer  system. 

This  will  enable  us  in  the  same  time  to  verify  the  results  obtained  with 
our  original  method  with  those  of  existing  programs. 

3.1.  Algebraical  anatysis  of  a  three-layer  system  (isotropic,  full  friction). 

3^  1  ^  1  ^_BoLn_cbr^__condltlons  _of^_  the_s^s  tern . 

We  consider  y*.,  ■=  p-a  a  °  ^ 

the  thickness  of  the  first  layer 


H2 

the  thickness  of  the  second  layer 

We  write  A,  for 

A,rnx 

b,  for 

b.wJ1 

C,  for 

C,  no 

1 ),  for 

and  \<  s 

E  *(•**>*  "A  Ea 

B7  (**)*«)  £-2 

Ei(^+)Al)  E3 

Boundary  conditions  at  the  surface  (z  =  0): 

r2.y>  ■. 

A,  +  £>,  =  * 

i 

Zr7  -O  • 

A,'  b,4-C(4  D,  ro 

Boundary  coitions  at  the  first  interface  (z  =  Hj,  x  =  mHj) : 

BaC'-*  +  a C2e-X x  P-j«rx 

Tf* :  A, t  -  B,e  \  C, 0+  *) ♦  5, (i-*) £' * »  A  jt  -  B C, (++*) e*  +  *  * 

W  A,t  -  +  <C,e  -x  3),  e.  *=  S>2 *- x -f  x  c2 e*  -  x  *3 

l*  ■  A  A  b.e*  V  c,  W[>  A  bi«:x  V  C2  0-*) *  *] 
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Boundary  conditions  at  the  second  interface  (z  =  +  Hg,  y  =  m.(H1  +  Hj)): 


(£•.  +  B>ae.y*  y  Cxty  *•  y 

rTX  ■  Alt’' -  £>2«ry  Ci  (A+y)ey  *  V>i (4 -y)t7  =  -  -  D-j 

W  •.  Aa*f-  batY  w;C iJ-yDa^  *  L  L"  -y  ^3ey] 

Aie7+E>2eV^C1(^y)eV_'D1^)t-'/a  L  [b  tT*- 

The  boundary  conditions  can  be  written  in  matrixform 
At  the  surface 

I  (  A,  B,  C,  I>,)T=  (A  o) 


where 


•C 


At  the  first  interface 

M.,  (A,&,C,I>,)T=  M,(AxB,C,Pa)T 

At  the  second  interface 

M3(Aie.,clp,)T= 

3ili2i_Sglut i on_of _ the_s^s tern_of _ 10_eguati  ons . 

We  start  from  the  conditions  at  the  second  interface  and  write 
(A1B,CiDx)'r=  Mj  M4  (ft,  D,y 
Matrix  M3  is  very  easy  to  invert: 
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that  we  write  as  follows 
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with 
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so  that 

-  4f  W"*.  *  ’V'X]  O  ^ 

We  notice  that  the  positive  exponent  ey  has  disappeared 

We  now  develop  those  matrixproducts  necessary  to  avoid  convergency  and 
overflowproblems.  For  the  others  we  only  take  notice  of  the  terms  equal 
to  zero. 

So  that 


Y\ 


w 


where  the  sign  +  denotes  some  constant  or  some  linear  function  of  m. 
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(A1TbxC1Dz)T=  -  %]f^.  -  (>h-L).Ul]  i^l>2>)r 


<?> 


The  term  (1  +  L).U^  will  be  part  of  the  constant  in  the  final  denominator, 
constant  whose  value  we  have  to  know  to  ensure  convergency  at  the  surface. 

We  now  consider  the  conditions  at  the  first  interface: 


(A,Ib,C,D,)Tc  M^1-  H2  (A2.li2.C2.  'Dz)T  ( 4 ) 
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We  again  develop  the  products  as  before: 
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because  it 
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verifies  that 

^U4l 
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^3241 

.  V\Ti. 

*o 

so  that  the  concerned  products  disappear  from  the  relation. 
One  verifies  also  that 

^3t4i  ^5Ml  =  * 

^3)42  .  Myi ;  2  -  Q 
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The  first  matrix  is  a  constant,  the  second  is  a  linear  function  of  y. 
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where 

W  =?  -  iL  (4~3Ly4^  ~  + 

^L,  ~  'Ll  (-4-^hV)J 
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w-.  Ajt*  -  3-O ^  - 

Lw  [A^e-  -  H>At  z  _Qi,  (2-4^4-*.)  t2- -  D*  (i-4 |»i,»-x)t7J 
u  •  +  V*  +  Cb  (a+x)  _  D3  (*-x)  c\ 

v  frACX  -V  CA  (*•**)  <-*  -  D4  ("i'O4  *J 

Boundary  conditions  at  the  bottom  (z  =  +  H2  +  +  H^): 

w.  A40  -  &ACk  -  C4  O^-O^-TU  (l-V^*  » 

T-^lPoayhc  •.  Ci,  =  o. 

3^2_.  2i_Sol  ut  i  on  _of^_  the_s^s  tem_  of  _  16  _egua  t  i  ons  ^ 

In  the  equations  of  the  conditions  at  the  third  interface,  Is  replaced 
by  its  value  taken  from  the  fixed  bottom  ndition: 

a  bAefc  v  TU  (1-4  +t)tt 
We  write  the  conditions  at  the  third  interface  in  matrix  form 

(_A»*>»C 4D»)t.  t^»o’T 


so  that 
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Boundary  conditions  at  the  first  interface  (z  =  H^): 

rz  A.c.%  b,£*  -  C,  +  D,  (^-2^, *■*)«:*= 

Xfz  :  A,e*  -b.e.'*  +.  C,  (2^,4-x)  «,*  +  D,  -*')«."*= 

Aa«-  -  1b  i.  4-  Co. c*  +•  1 7  ^ 

^  •  A,®.*-  13,<TX  -  C,  (2- 4)*,*-*)  a*  -  T>, 

Fv  [A2 ^  -  £>a*- x —Ci (a-2^T-x)t*-  t>2  (2-^i 4-x)c*J 

Ia-  A,4-X  +■  ^),C  4-  C,  ^-44-X^  eJ*  —  'D, 

^"lA.  4-  U>2e  +  —  ^2  C^~x^  ®"  J 

Boundary  conditions  at  the  second  interface  (z  =  Hj  +  H^): 

(Tx:  AjO'  +  £>247  -  C-t  (^-2^1-7)^  v  D2(>i-2vv^ y)*-  y= 

+  b34^  -  C-2,  (/1-2^,,-y)  «-*  4  ^  ^ 

trx:  A-zt^  -  bae7  ^-01(2^4^^  4  ‘I>tCAW*“)Oc- '^=‘ 

A»^  -  4-  C  3  ^2y>.  v  1)3 ^ 

w:  Aj*y  -  biCy  -  Ca  [x-Ayi-y)*^  -  %[p>-4yi  4y)e  ^  = 

Vw^A^e.^-  ^3®^  -  C -3  -  5>3  (J2-4yj4-'0®',J 

U;  A i^y  V  b^€^'!|,  -4-  C2 

Ww  ^Aje?  V  bse^  4-  Cy  C^+'l)  ^ 

Boundary  conditions  at  the  third  interface  (z  =  Hj  +  +  Hg): 

Tz  ■■  Ajef  +  ba^i  -  C >  [a- 2y  3-*)^  4.  X13  [/<-2y3  ■>‘x)e.  - 

As^*  v  -0/,^'2^.4-a)  t2  4-  3U  (•4-2)Ai,4-z')  £ 

t-ta  ;  A*e*  -  b^eT*  4  C3  (2^3  +  20e-  4  = 

Anc*-  f  C*(2n*x>*  4  V*  (j2.^-*)eX 
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3.2.  Algebraical  analysis  of  an  isotropic  four-layer  system  with 
fixed  bottom. 

3i2ili_Boundar^__condi tions^ 


In  this  analysis,  we  choose  as  fixed  bottom  condition  the  one  described 
in  §  1.2.2.  (w  =  0). 


We  write 

A,  -  Ajwi1-  £>,  =  C i  ^ 


■s  2)-,WJ 


Fw«  Eil*J±±l 

E  2  L>,+  Vw') 

Lw  *  E> 

E  h  C>,+V"‘0 


ru  =  ">  - 

ww  .  w 

E-j 

Lu_  .  v,  M44^ 

E4 


X  -  w*V\, 

y  »  v 

y  ,  W>  [H,  V  Vi*  +  fU) 

U  ~  -2t  f  2 


where  Hj,  H^,  H3  and  are  the  thicknesses  of  the  four  layers. 


Boundary  nditions  at  the  surface  (z  =  0): 

A,+  Bv 

Xrx  5°  ^  t  ~  B  *  4-  O-t  •  2.  y*. ,  2^,  a  Q 


P-  1.00 


El 19 100000  U1=0.D0  H 1  =  1 0 . 00  A=  10.00 

E2-  10000  U2=0.50  H2=20 . 00  R=  0.00 

E3~  1000  U3=0.50 

*********************************  ******  a  ***  *********  ****************  a*  *  a  a  a  «  a  a  a  a  a 

*  a  SURFACE  *  BASE  *  SURFACE  *  BASE  a  SURFACE  a 

*  *  COUCHE  1  *  COUCHE  2  *  COUCHE  2  a  COUCHE  3  * 

**********  *  ***************  a  aaaaa 

•  CONTRA I NTES  VERTICALES  *  *  +0.247002*  +0.247002*  +0.029472*  +0.029472* 

*******  *********************** *********  a  a  ************ ***************  a  »a****a  a  aaa 

♦  CONTRA I NT  ES  RADI ALES  *  *  -1.B64850*  +0.035817*  -0.204340*  +0.006091* 

*  a  *******************************************  *****  *******  a*  *************  a  *  a  a a  « a  * 

♦CONTRAINTES  CIRCONFER.  *  *  -1 . S64850*  *  -0.2043+0*  +0.006091* 

*********************************************************************  a  ********** 

♦FLECHES  *42. 183E-03*  *  *  *  * 

a******************************************************************  a  *******  **  a  a  a 


Table  1.  Stresses  in  a  three-layer  system  (BASIC  Program). 


RESULTATS  DANS  L  'AXE  DE  LA  CHANGE. 

El*100000.  UlsO.50  MlslO.OO  A=  10.00  p*  1.00 

E2*  10000.  U2»0 .50  H2s20.00  R*  .00 

E3*  1000,  U3»0.50 


SURFACE  BASE  SURFACE  BASE  SURFACE 

CCUCHEl  C0UCHE2  C0UCHE2  C01CHE3 


CONTRAINTES  VERTICALES  *  .207002  .207002  .029072 
CONTRAINTES  RADI  ALES  *  -1.869850  .035B17  -.209390 
CONTRAINTES  CIRCONF.  *  -1.869850  .035817  -,2093«0 
FLECHES  .0022  *  *  * 


.029972 

.006091 

.006091 

ft 


1 


Table  2.  Stresses  in  a  three-layer  system  (FORTRAN  77  Program) 
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3^1 i4i_Comgarlson_w1 th_exi§ ting_grograms . 

First  we  have  computed  the  stresses  in  a  three-layer  Isotropic  system  with 
a  program  written  in  BASIC  and  for  which  the  equations  had  been  developped 
in  complete  closeform  (VAN  CAUWELAERT,  1983). 

The  accuracy  of  this  program  had  been  checked  earlier  with  the  results 
published  by  JONES  (1962  ). 

Then  we  have  written  a  new  program  in  FORTRAN  77  based  on  the  developments 
presented  in  paragraphs  3.1.1.  and  3.1.2.  and  calculated  with  this  program 
the  stresses  in  the  same  three-layer  system. 

The  results  obtained  by  the  program  in  BASIC  are  given  in  table  1,  those 
obtained  by  the  program  in  FORTRAN  77  are  given  in  table  2:  we  notice  that 
the  agreement  between  both  is  perfect. 

The  complete  listing  of  the  FORTRAN  77  program  is  given  in  appendix,; 

The  matrices  are  dimensionned  and  upbuilt  in  the  beginning  of  the 
program.  The  variables  x  (=mHj)  and  y  (=mHj  +  mH^)  are  introduced  in  the 
matrices  by  instruction  214.  The  products  between  matrices  are  carried  out 
in  the  same  instruction  214  by  calling  subroutines  800,  810,  830  and  860. 


I.R.  41. 


§;i;3i_8®l§9i2!]_?2r_£!3§..¥§rti£§!_dffl§ct1on_at_the_surface. 

We  notice  that  the  values  of  the  parameters  given  by  (10)  and  (11)  all  converge 
normally,  except  for  parameters  and  which  contain  each  a  constant  in 
the  numerator: 


4.  L)  B4  »  %  (**W)(-«L) 


4(-hW}(/ivL)  +  k&ix 


L^V<) L)}  •  eti2 

V  *  V 


The  influence  of  this  constant  can  be  eliminated  as  indicated  in  §  2.3. 
Here  it  is  nevertheless  easier  to  eliminate  B^. 

From  the  boundary  conditions  at  the  surface,  we  have 

=  a  -  A  ^ 

The  relation  for  the  deflection  at  the  surface  is  given  by 


/*Q 

w.  -b«x  *_thL 

4  m  E, 

lA«  ■  ^’3 

>»+u. 
=  £7 

I  \T, 

)  M 

'o 

[a  -  2A.) 

6 

1 

ll 

)v> 

KJ-\  A.cU*  0^ 

•i  K  J 


T 


E, 


Vvi 


The  first  integral  of  (12)  is  solved  analytically,  the  second  one  converges 
fast  and  safely  during  the  numerical  integration  procedure. 
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The  parameters  Aj,  B^,  and  0^  are  then  obtained  from  (6) 


A  ~  -  ± 


B'  =  x 


C  =.  -  d. 
I  A 


(>i+V <)e-  (V)  ^3  +  O’1) 

>^>i+L)c  jj’liiaiL  O'*)  +■  ^iiiiL  O'1) 

-2{y-x) 


•’I 

0*V)O+L)  +  «-  ^  ^[M,24i(*.»)&3  +  Mi2AiO*2)^]j 

(A+*)lly  [tt**,  (v )  fc3  4-  (3.a) 

v^vL’)^  [miu2L(3,0  e>3  +  MnizL  (*,*}  DaT]^ 

A  =  ^*|^+W)04-L)P3  t-  ^jV|az„  +  Hi4t 


Parameters  A^,  Bgj.Cg  and  D2  are  obtained  from  (3) 

A2*  -  J^V  &3  ♦  H3^1 

£>1  *  X  C'1* 

Cx~*‘>  31  ill  •  +  M3141  C3*1)  D3J 

3>i=  |(^^Tb 
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I.  M- 


X  •  M  i-2^» 


X.  miL  “ 


The  terms  a..,  b . .  and  c.  .  are  the  results  of  very  simple  computerprocedures. 

'  J  *  J  *  J 

We  can  now  write  (8)  in  extended  form: 

_ly  _2x 

I M+k)  (4*U)  -  (_>»♦  V)  .a„  e  4.  b,,.*'  -  [a+L).Cu  «-  J  t>3 

.  -lv  x  -^Y-*)  .  . v  -2*“\  „ 


_2V  -2(y-xJ  j.2y"\  n 

j-  ^+k)(-i+L)  -  ^+W).  V  *  J  K3 

-(A4-k).ft2i ^  -^+L).  C22«-  D3 


We  write  this 


[(^WX'1+L) +  +  d,a  t>3*  ^ 

^^k)(-l+L)+dx,]E>3  +  +  *  ° 

where  (1  +  K).(l  +  L)  is,  as  a  matter  of  fact,  a  constant  and  all  d^,  d^, 
d£^  and  d^  converge  normally. 

Finally  one  obtains 

D  4(a+)*)(*+V)  *  AA12  ^  4  -  4  d 

Di=  — * - - -  i'jo  - - 

V  V 


where 


y3  4.  ^4.Vr)(-n+9^«  ♦<*iz  *<*«)  *  d*  <kzi-  cAu <*ij 
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•  ^314 1 


O  o 


+  + 


o  o 


^1122.-  ^L 


4-  + 

0  o 


4-  4* 


o  o 


-  n, 


—  ^>H2L 


-  O 


so  that 


(7\,6,c,  Di)T.  A- [  (-i+V)(a+L).  Ui.  -  vM3iai 

+  -(^L)e‘2'M,lllLj(B3  ,3>T 

The  positive  exponent  e2x  has  again  disappeared  together  with  the  highest 
negative  exponent  e~2^y  +  x). 

Finally  we  consider  the  conditions  at  the  surface 

X(A,e»,C,5,)T  -  (A  0)T 

We  replace  in  (7),  (Aj  Bj  Cj  Dj)1  by  its  value  from  (6) 

[_2y  -2ty-x) 

-(^i+L^e  MjiiilJ  [&2>  =■ 

fH.W)(i+L)X.Uu=  f-.+  uj^+L)  (  ^ 


/ 
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We  replace  in  (5),  (Ag  B2  C£  D2)T  by  its  value  from  (3) 


(A,b,C.D,)\  ±.  w).Uw  +  +  e.  H127l]  • 

\eIyM„V.  -  C^L>UlJ(E,5D3)T 

«  Uk  Ul  -  Uk  Mj,/,! 

-'zfytJ')  -2[y-*X  *>* 

-V  e.  •  +■ e  Mvjxi 

-(hL'^JWUl  J>0T 
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The  product  we  are  interested  in  for  later  computation  is 


This  product  contains  a  constant,  a  linear  function  of  y  and  a  linear 
function  of  z.  The  term,  function  of  y.z,  a  nonlinaer  function, disappears. 

The  other  products  are 


\ 


I 


l 
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Next  we  write  the  conditions  at  the  first  interface  in  matrixform 

■■  j  ■ 

(A,  1)>,  C,  D,)T  =  Vi’’.  (OtKCtJi) 


I.R.54. 


The  first  matrix  is  a  constant,  the  second  is  a  linear  function  of  x. 


,  .  .  M 
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•V 
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o  0 

O  +• 

O  0  / 

o  o 
+-  o 


=  M 
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o  o 
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-  n 


U2I 


/ 


A  - - -  \  TV  Mi.  ► 


-1 


_2y 

t  c  .Ni,  v  t  .Vu 


-Xtt-y) 

+•  Nil 


-2  (*'7') 

vt  1  Vv^  +  e.  •  JV23 


-XV  Va*  v  *  2*  Wi4 


f  e  "*?  Ni2  +  1  l  $4  ^0 


T 


i 

*■  Mail  *- 
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The  products  which  will  not  converge  normally  are 

elx  «.  ^  .  Mtu» 

4T^*'y)  Han  M-jui  tijui 

«x  Mni.  (Vln-M*  + 

because  that  it  can  happen  that  x  \“-y 

or  x  >  z.-y 

or  x  f-z 

and 

M.H,  (m«  m‘0 

One  easily  verifies  that 


Mat!  .Mu)  -  M?iti  «« 

M\ixi  •  (_^v  v  Mao)  •  Mya-  Mu.  *  *» 


so  that  again  the  concerned  products  disappear  from  the  relation. 
The  product,  we  are  Interested  in  for  later  computation  is 


|Vln  Mi,  *  nv^xJJ^a.vU,  - 


V(.K,-Wr)| 


t>  o  \ 
P 


(.LrU) 


O  o 

©  o 

0  o 

o  F,  W, 

O  O 

O  © 


/ 


=  T. 
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FVWVLS  +•  FvV-il-i  +•  ^-4^J  -  V-iW^S.  [a- 


One  has 
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One  verifies  that  the  products 


^1111'  ^  L  ”  ° 

TV  \  « 


The  other  products  are 

^>r»  — 


and  finally 


(p>,Tb,C^0T=:  “ 


( >4  |a-  0 


r  -i-i  -ly  _  -2-1  T-y)  -1l2->')  X, 

l  e,  *,  .  Ta,  +•  *  •  Txx  +  T*  r  *- 

-2-V -2x_  -Z(^V-t)  -2*  -X(V--y)  <— 

VI  •  '1^  *  «-  '2.4  +  *■  >14  *  e  *-  •  *»2. 

->  -7*  -z[y-zU  .  TL 

-2 ly-*)T  T, 

v  «-  «.  . T^,  +  «-  .  >42  ■»  e  •  '2*3 


v  ^  T/,4  -r  T«  J[t>4 


(fnl.3) 
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Finally  we  considerthe  boundary  conditio  s  at  the  surface: 


(4-V)  (H- V)  \  T  ,  .T 

ly,  j 


that  we  can  transform  into 


64(A"h'XvH(>l'h) 


»  -  O'  VO  CA'V0 


\  -l 


c  X,  L\1 


£-*t 


Oil  X, 


*«  1  o„  mt=  (»  «y 


Extending  matrix  (d^),  equal  to  Tu,  we  obtain 


G^t  (■*-  {* 


! _  (  '  ' 

0  \  ' 


-  v)  O'  v«)\[  CZ«  CZl 


C»  C,i 


Cil  Cyz. 


Ca,  cW 


/  F.V.L,  ^,W|L,CA-A>^  ►  \ 

\  -T.k,  L,  \W*  -  nVfxU]  / 

+-Ix(>rl-?i).WvW  »-y(k,-wv)F,u  *  *(l,-u)f,w,]  (  0  \  iUh)Ts  0 

1  ) 


This  last  relation  can  be  writtenin  extended  form 


U„  *■  4-  (>i- 2^,)  C.*,  4-  F,  '-r.L-.J  $4 

v  Cu_  —  1^.,)  Cjj_  4-  2^*i)  Ci^i  +-  -L  ^F,Vil-,  *- 

u[Frh)WU  v  v(W.-Wt)T=,W  ♦  z(L,-li) F,W,  ^  g  ^‘>^1 


Cl,  -  c  lt  V 


2^_,  c*,  v  2^.4  c^,  -  TiW.L,^  B/, 


*■  |  c,r  -  clx  v  2^,  c-^x  4-  2^,  c^x  -  ]jF,Vc,  L,  -TzV-»  Lz^ 

->(F,-r^  Vr^Li  -  v(Vf>-V<x)FvL1  -  z.  (_L*-Lx)  F,V,^  =1  <a 


We  can  write  this  system  as  follows 


^rvj|  4-  FtK  L,')  ^  V\]X  4-  j^Fk\L.,  *" F^Vvl-iJ  v  ^(xyl)|*04  3 
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(kx,  -F,W,  L,')  Ba  *  j  ^xx  -  ^F,V\L,  (4- 4 )*j,)  TrWW^  -f 

where  all  n^j  converge  absolutely  normally  and  f(xyz)  is  a  linear  function 
of  x,  y  and  z. 

Solving  this  system,  we  obtain 

Bii=  _  L  [?,w,L, (^-Ay-i.) -T.  WU}  -fin1)}-  '? 


I  - R .  59. 

where 

V,  O'*  ♦  >**•) 

_  *M  1.  ^V.u,  [>*-  A^i,)  -TxWtLi^  ..  viii  4p  (>'-^hA)  v  Tt  WU] 

-  *,Wr,U(  - 1-  j>,  V,L  (>i-  Am)  -TtWJUj  *-  T.W.U  ^  j^V.L,  h'^K)  V^WtL^J 

4-  F.W,1-.  f-(r*x)  -  T.W,t, 

v-  ^  (*„,**.>  v'1,>  ***)  v  Kv,M  ^xU 

The  linear  function  f(xyz)  has  disappeared  In  the  denominator. 

The  function  f(n-jj)  converges  normally  so  that  the  limit  value  for  the 
denominator  is  a  constant: 

JL\aa~  V  «  T\  V»  L,  ■  W  Lx 

m  — >  *o 


I 


3.2.3.  Values  of  the  parameters  A..,  D... 


At  the  bottom  of  the  first  layer,  parameters  A^  and  are  factors  of  the 
positive  exponent  ex.  Parameters  Bj  and  D1  are  factors  of  the  negative 
exponent  e'x. 

Thus  when  computing  stresses  in  the  first  layer  we  have  to  input  at  least 
two  positive  exponents  which  necessarly  will  lead  to  overflowproblems. 
Therefore  we  express  in  the  program  next  modified  values  of  the  parameters 
Aj,  ,  Cl  and  A1.ex,  B^,  C1>ex,  D^. 

Then  if  we  have  to  compute  a  stress,  let  us  say  at  a  depth  2x/3,  we  multiply 

x  x  -x/3 

the  parameters  A, .e  and  C,.e  by  the  negative  exponent  e  and  the 

1  1  -2x/3 

parameters  Bj  and  D1  by  the  negative  exponent  e  '  and  avoid,  in  doing  so, 
all  overflow  problems. 


The  values  of  the  parameters  Ap  B^,  and  are  given  by  (sol.  3)  of  §3.2.2. 
To  obtain  the  values  for  A^.ex,  B^,  C^e54  and  D1  we  split  (sol.  3)  into 
two  parts: 


(*' 


o  C,  o) 
A 


(  o  B, 

L 


o 

J  c*  »oT 


We  notice  now  that  the  matrices  T^  and  T^  contain  nothing  but  zeros  in 
their  first  and  third  rows,  so  that 


and  that  the  matrices  T^  and  T^  contain  nothing  but  zeros  in  their  second 
and  fourth  rows,  so  that 


o 
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Thus  the  values  of  and  depend  only  on  the  matrices  and  . 
This  has  to  be  so,  because  only  the  exponents  multiplying  the  matrices 
and  T^  do  not  overflow  when  they  are  multiplied  by  ex. 

The  exponents  multiplying  the  matrices  T^.  and  T. .  could  overflow  when 
multiplied  by  ex:  for  example,  the  exponent  e"^  ~  multiplying  the 

matrix  T 14  if  2(t  -  z). 

We  obtain  then 


ex  (A,  o  c,  *)  -  _ 


■(2y~x)  ~a(_r-z.) 


A  (V {a-y-l) 


1  V.  <  -v_  -J  -1-7  )  - I  -J  __  -T- 

•  «■  .  T-,  +  e  ■  t;„  -v  e.  ,t;x  +  c.  T, 


-{ir-y) 


A-  -y-  -x  -‘l(*-y)-r  -*  -y 

+  e  e  U3-we,.ev  .  T, 


2-i  T  •  '22 

-y  -z(v-z). 


2  A 


'3 A 


A-  «- 


£>,  o  3),^  _ _ 

U-^) 

r  -z(?-7)  '2[h-z)  -2[y-y)  -2[h-z) 

Le  ^2  ♦  T,  *  «-  Ti,  ■'* 


+  «■  TA2  a  e  -  >a4 


®o' 


We  follow  the  same  procedure  for  the  param  lters  A 2»  Bg,  C2  and  D2  and  obtain 
from  (sol.  2): 
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2) 


■V 


•'tt-y) 

«■ 


,  -K  \  T  ^ 

C°  0  *v  ’ 

ft  N„  >  »  *’l>  "V*  ♦*..](>*  C*)r 
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Finally  we  obtain  for  the  parameters  A^,  B^,  and  D^sol .  1): 


2.  /  N '  J\ 


r«^  •  ^riti  *■  *■  [fa 


/  O  O  $3^  s  - - .  t  .  VJ)5l4( 

^-v~0  L 


( bi  »*) 


3 . 2^,4.  _The_deflect  ion_at_  the_  surf  ace . 

The  vertical  deflection  at  the  surface  is  given  by 

w=  1  -^4  1±i  j^A , b-v)  - 23)' 

The  parameters  Bj  and  do  not  converge  normally  because  they  contain  a 
constant  in  their  numerator. 

Utilizing  the  boundary  conditions  at  the  surface,  we  express  then  ^  and  D1 
in  function  of  Aj  and  C^: 

a  />*  2A(  t  C;  (-4-4)^ 


so  that 


A,  -  2-C,  (^-2^,)  -iD,  O'1-2}**)  ” 

\^a~  2A,  ac, 

w-  f  J»0"Q-TM 

e  ,  /.  *-> 

_  <K).t»r3.|»r).TH  j- ^ 


The  first  integral  can  be  solved  analytically  (§  2.3).  The  second  one 
converges  absolutely  normally. 
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2i?i§i_I!]§..§£!;§s§§§_§D<L^i§B!?S§roents_in_the_first_layer. 

We  compute,  for  example,  the  vertical  stress  at  a  depth  0<  H  H^, 
Its  value  is  given  by 


I  3,  (*,►)  A,*-"  V  ft, 

/»  u 


_^>i 

e. 


-w.H- 


-  Cr«. £>»-  2^,-u^H) .  e.  v  J  eU^ 

If  the  value  of  H  is  very  small  (if  we  compute  a  stress  near  the  surface), 
the  values  of  B^e"111^  and  D^.e'mH  will  converge  very  slowly. 

We  know  from  the  preceeding  paragraph  that 

B,e.  a  A,  +  A)-.  C,  (>i- l|*v) 

5,Cw,M  -  SA1«-'~H  ♦  c, 

So  that 

(4+u.^ft.  _  2A, a-  C,  (*- 1^.  vkH 

We  write  then 

X  M.  % M )  A,t. t  C.J (*- 2,.  ct-x-t*h) 

i.  L 

*(nvJ0cta,H- 

We  can  again  split  this  integral  into 

\A  u  H 

CTa,»  1  ^po.  T*a  Cyy,r)  Cvy,°')  (/»  +  l*»H)  C  tU^ 
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i  •  •' 
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The  first  intergal  can  be  solved  analytically  (§  2.4),  the  second  one 
converges  normally. 

Of  course  the  same  procedure  can  be  applied  for  the  computation  of  all 
other  stresses  or  displacements , in  the  first  layer. 

3^  2_.  _Jhe_  s  tre  §  ses  _  and_dl  sg]  acement  §  _  In  _  the_  second  _l§yer . 


Normally  one  should  not  have  numerical  problems  in  the  computation  of 
stresses  and  displacements  in  the  second  layer,  except  for  the  case  that 
the  first  layer  should  be  very  thin  (which  can  happen  with  overlays). 

The  terms  B2;e”mH  and  D2.e_rTlH  could  then  again  converge  quite  slowly. 

The  solution  is  then  obtained  as  follows. 

We  write 

d  .  B  - - ■  •  Mm  (z.i  ,  2.z)(S>* 

_Da-  . 

where  N^(2.1,2.2)  are  the  first  and  the  second  term  of  the  second  row 
of  matrix  N^. 

Dx,  -  Da - - - —  .  Maa  ®a)T 


bAy  =  B/,  -  fr-HB-Hfr-M)  S  -f ^ 

b  »  ^ 

-  -  B*z 

.  t.v.a. 
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-  D/,  v 
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T 
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we  split  then  (B^  D^)T  en  write 


e.,=  b„  *  — ^-77 — ->v.(n,i->)k  5«)t  ♦  (V  C«)TJ 

•y 

=  bi,  *  baa  +  - - - :  2.2)  (^2  5>4i) 


Da  -  Do,  +• 


— - Pi>0T  ^  ^Ojl 


=  Da,  V  Daa  +  - - - Hu  (4.»,4x)C^i  P*0 

We  notice  that  the  terms  82^,  Bgg*  and  D22  converge  normally. 


!V„  (21  ,  22)  (  ^A2  D«)  s. 


\*-'  d_  j"  T,  ^L.  (^-4^)  -  T?  V^Lj^  -f  -CX7Z^  y 
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=  W,L,  V-,L2  ^  4  X  (T.-^  Y 


l)a-x.=  - - - -  ton  (4. >,42)  (Titoi  DaO  =■ 

K(a-^)(^) 


»  ±tl.  V<aLaT,V,  L, 

4  V 


We  split  now 


7=  7I+  b,W,L,.T2WTL: 

where  V,  Converges  normally- 


We  divide  then  the  numerator  by  the  denominator  and  obtain 


0  t  b-  *  wO  -^0  +  xt'F>_?0 


P,<W 


The  first  term  can  be  integrated  analytically  (§  2.4),  the  second 
one  converges  normally. 


?;?iZi_T!]§_§Tr;§sses_and_disglacements_in_the_third_layer. 

Following  the  same  procedure  as  in  the  preceeding  paragraph,  we  obtain 


^  O-k)  V| 

^  W-.F,.  K-xFj  .  v 


-'U 


References . 


8RUMISTER  (1943):  The  Theory  of  Stresses  and  Displacements  in  Layered 
Systems.  Proc.  Highway  Research  Board.  23  rd  Meeting.  Washington. 

JONES  (1952)-  Tables  of  stresses  in  Three-layer  Elastic  system. 

Annuel  Meeting  of  Highway  Research  Board,  1962.  Washington. 

Me  CULLOUGH  (1981):  CRCP  Design  based  on  Theoretical  and  Field  Performance. 
Proc.  of  the  2nd  Int.  Conf.  on  Concrete  Pavement  Design.  Purdue  University. 
INDIANA. 

TIMOSHENKO  (1970):  Theory  of  Elasticity.  Me  Graw  Hill  Book  Company.  TOKYO. 

VAN  CAUWELAERT  (1983):  L'fclasticite  anisotrope  appliqu§e  a  la  mScanique  des 
milieux  granulaires  et  des  roches.  PhD  Thesis.  Ecole  Polytechnique  FSdSrale 
de  Lausanne.  Switzerland. 

WATSON  (1966):  A  Treatise  on  the  Theory  of  Bessel  Functions.  Cambridge 
University  Press.  Cambridge.  Great  Brittain. 

WAYLAND  (1970):  Complex  Variables  applied  in  Science  and  Engineering. 

Van  Nostrand  Reinhold  Company.  New  York. 

WESTERGAARD  (1926):  Stresses  in  concrete  Pavements  computed  by  theoretical 
Analysis.  Public  Roads.  Vol.  7  nr  2. 


APPENDIX 


Three-layer  elastic  system 


Computer  program  F77 


UR  PART  IE  :  DECLARATION  DES  VARIABLES 
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Appendi x/1 


FUHMATC+PRESSICN  APPUQUEE  EN  SURFACE 
READ (5# '(F5.2)*)F 


M(l)*-2.«0ETAl*Xl/t'/El  (0) 
h(2)*UETAl+CETBl*XltDETCl*X*CETD!*X*Xl 

M(3)«-.5*CDET»l+DETBl*Xl+{3.AX)*DE1Cl-(3.-X)*DETDl»Xl) 
H(X  jsOEt  A2*Xt ♦DETB2*X1 ♦DETC2*X»  X6  +  DETD2*X*X 1 


6EHE  P»RT1E  t  IPPPESSIO*  DES  CUMHAJMES  ET  OES  DEPLACEPENTS 
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